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Abstract 

We argue that precise measurements of charge and magnetic radii can meaningfully constrain 
diquark models of the nucleon. We construct properly symmetrized, nonrelativistic three-quark 
wave functions that interpolate between the limits of a pointlike diquark pair and no diquark 
correlation. We find that good fits to the data can be obtained for a wide range of diquark sizes, 
provided that the diquark wave functions are close to those that reduce to a purely scalar state in 
the pointlike limit. A modest improvement in the experimental uncertainties will render a fit to 
the charge radii a more telling diagnostic for the presence of spatially correlated quark pairs within 
the nucleon. 
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I. INTRODUCTION 


Are there correlations among qnarks?” is a simple yet central qnestion in the s^nd 
of hadron strnctnre. One conjectnre that has received attention over many years 


s, 


and renewed advocacy in recent years 


Andy 


, is that pairs of qnarks form identihable 


“diqnark” snbstructures, where diqnarks have a smaller size scale than hadrons overall and 
have dehnite spin and isospin qnantnm nnmbers. The name “diqnark” has been used with 
various meanings. It sometimes refers simply to the two quarks remaining after one quark 


is scattered out of a baryon j^; of more relevance to the present investigation, it has been 
used to refer to a pointlike or nearly pointlike state with the quantum numbers of two 


quarks [1|, |3[ . We wish to learn if experimental evidence indicates the presence of diquarks 

that are small, but not necessarily pointlike, two-quark correlations. At some level, the 
existence of quark-quark correlations has been known for decades, since the nucleon-A(1232) 


Q 


mass splitting is generally attributed to a spin-spin interaction [Tj that is attractive for 
spin-zero isoscalar quark pairs and repulsive for spin-one isovector pairs. The attractive 
interaction shrinks the size of the favored quark pairs, but is unlikely to force them to be 
near the pointlike diqnark limit. 

In this work we are interested in nucleon form factors, which provide direct information 
on hadron structure. Nucleon form factors have been studied previously in a diquark-quark 


model j8[, but the diqnark has in these studies been taken as a distinct state, sometimes 
including its own intrinsic form factor. Other recent work on the consequences of diquarks 


in hadron phenomenology can be found in Ref. P|. Here we start with a three-quark state 
and describe the nucleon with wave functions that contain quark-quark correlated states 
bound to a third quark. We include scalar diquarks, which refer to isospin-zero states, 
which must also be spin-zero if bound in a relative S-state, and vector diquarks, which 
are the corresponding spin-one isospin-one possibilities. We do not overlook the role of 
Fermi statistics. If the diqnark state is very small, the effects of antisymmetrizing the 
state with respect to interchanges of one of the quarks in the diqnark and the third quark 
may be neglected. Many authors have yielded to the temptation of assuming a simplihed 
hadronic state that is not totally antisymmetric. However, if the diqnark is not much smaller 
than the hadron, overall antisymmetrization is important. All of our states are completely 
antisymmetric among all three quarks. We focus in this note on the charge and magnetic 
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radii of the proton and neutron. These are low-momentum transfer quantities and allow us 
to use nonrelativistic models. Further, we limit ourselves to three-quark Fock components of 
the nucleons, without relative orbital angular momentum. Our wave functions are written in 
Jacobi coordinates, which naturally allow a quark-quark correlated state of one length scale 
bound to the third quark with a wave function of a different and presumably larger length 
scale. Scalar and vector diquarks in an arbitrary mixture are allowed. The formal details of 
our calculation are given in Sec. im and the results of our hts, including the allowed scalar- 
vector diquark mixture and the ratio of length scales, are presented in Sec. IIIII Section IIVI 
summarizes our conclusions. 


II. FORMALISM 


The nucleon wave function is completely symmetric under spinxflavorx space, with color, 
as usual, carrying the antisymmetry necessary to guarantee the proper fermionic statistics. 


The flavor and spin wave functions are of the standard textbook forms jUl, for example. 


la. 


Pms = -^[^uud - {ud du)u] , pma = - du)u , 

xms(T) = -^[2 in -(Ti + U) T] , Xma(T) = - U) t • 

The mixed-symmetric (MS) wave functions are isovector (vector) in flavor (spin) and are 
symmetric in the hrst two quark indices. Likewise, the mixed-antisymmetric (MA) wave 
functions are isoscalar (scalar) and antisymmetric in the hrst two quark indices. One may 
construct a nucleon state that contains a pure scalar (vector) diquark by multiplying PmaXma 
(PmsXms) with a wave function 0 ms— either in position or momentum space—that is also 
symmetric under exchange of the hrst two quarks. Such a combined wave function is, 
however, not completely symmetric under all quark exchanges because it treats the hrst two 
quarks as special; one must add to it a similar expression that distinguishes quark pair {13} 
and one that distinguishes quark pair (23). 

We denote the completely symmetric wave function obtained in this way frompMAXMA^MS 
as IF), and the symmetrized wave function obtained from PmsXms^ms as \V). Signihcantly, 
in general (F|ld)7^0. Taking |1) to be the normalized form of |F), we dehne a state |2), via 
the Gram-Schmidt process, that is proportional to \V) minus its projection along \S), so 
that (1|2) =0 and (1|1) = (2|2) = 1 . Finally, the full wave function 10) consists of a linear 
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combination of these states: 


1'^) = cos6* |1) + sin6* |2) . 


( 2 . 2 ) 


An overall sign flip in this wave function is unobservable, and therefore all observables 
obtained depend only upon quadratic forms in cos 6*, sin6', or equivalently, upon sin 26* and 
cos 26*. 

We employ a notation that can be applied to any choice of the nonrelativistic 3-quark 
spatial or momentum space wave function. We begin by assuming only that the momentum- 
space wave function in the center of mass frame may be written as a function / of the form 


4,(abc )= /(pyy . 


(2.3) 


where {abc) represents the three single-particle coordinates in any hxed order, and the mo¬ 
menta in Jacobi coordinates are dehned as 

+ Pb - 2pc) ■ (2.4) 

Note that (^{abc) =(f){bac). 

First we dehne the basic overlap integrals 


M = {(f){abc) I (l){abc)) , 

p = {(j){abc) I (j){acb)) . (2.5) 


Of course, if 0 is normalized then M = l. Moreover, due to the symmetry under coordinate 
exchanges, every other such overlap integral is equal to M, p, or p*. 

We are interested in matrix elements 


(0(a6c) I Od I 0(e/^)) , (2.6) 

where in conhguration space, Od = exp{—i Q- r^). It is straightforward to show, again owing 
to symmetry under coordinate exchanges, that only four independent matrix elements arise: 

Ai = { 4 >{abc) I Oa I 4 >{abc)) , 

A 2 = { 4 >{abc) I Oc I 4 >{abc)) , 

61 = {(j){abc) I Oa I (j){acb)) , 

62 = {(j){abc) I Ob I (t){acb)) . (2.7) 
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Clearly, Ai^ 2 (Q = 0) = M and 5i,2(Q = 0) =/i. All other permutations of indices lead to matrix 
elements that equal either Ai, A 2 , (5i, 62 , or 62 - In the case that (j) is manifestly real, 

so are /i and 5 i^ 2 'i otherwise, the following expressions remain true under the replacements 
—>• Re/i and 61^2 Re(5i_2- 

Using Eq. (Q, these four integrals may be written as special cases of the single master 
integral 

X /(l(cos7)pL“'’''^+ (sin7)pf''‘^^-y^(sin5)Q|, |(sin7)pL“'’''^- (cos7)pi''’'^^-y§(cos5)Q|), 

( 2 . 8 ) 


where Ai=X(0,f), A2=X(0,0), 5i=X(-f,-f), and 62 =I{f, 0 ). 
The nucleon electromagnetic form factors 

Ge{Q^) = , 

i 

i 

where the Cj are quark charges, can be expressed as 

fisinO \ Ml 


G{Q^ 


= - cos U - -p _ , 

3 2 ^{M + 2 fi){M - fi) J 2M + /i 

fism9 


sin 6 


— cos 9 - , _ , , _ 

Sy 2 ^(M + 2/^)(M - 11 )) 7 (M + 2fi)(M - fi) 

{2M + fi) sin^ B 


M. 


36(M + 2/i)(M — nY 

ds , 

where the matrix elements M are given by 


Mip; = 2Ai + 4A2 — + A 62 , 

Mie = 2(Ai — A2 + (5i — 62 ) , 

■^ 2 ,E ~ 3(<^i + 252) , 

MIe = 0, 

■^ 3 ,E ~ 9(2Ai + (5i) , 

M^e ~ ~6(Ai — A2 + — ^2) , 

■^l,M ~ '^2^2 + + 4(52 5 

Mi m ~ ~2(A2 + (5i + 62 ) , 

MI 2 M ~ 2A1 + (5i + 662 , 

MIm = -2(Ai + 2(52) , 

-^3,M ~ 3(4Ai + 5(5i) , 

■^ 3 ,M ~ —2(2Ai + A2 + 5^1 + 62 ) 

From these forms one can show that G^^{0) = 

= 1 and Ge(0) = 0, independent of 9 

radii are dehned by 

2 6 dG{QY 

^ - G(0) dQ2 > 


(2,9) 


( 2 . 10 ) 


( 2 . 11 ) 


charge 

( 2 . 12 ) 
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with the denominator G(0) omitted for the electric form factor of neutral particles. 


III. MODELS AND RESULTS 

In the models used in this paper, the scales of the momenta p_ and p+ are set by param¬ 
eters a and (3 respectively, each having dimensions of length. All dimensionless observables, 
such as the ratios of charge radii, are therefore functions only of the two dimensionless pa¬ 
rameters a/j3 and cos6'. The small-diquark limit is dehned by a//?-Cl and the no-diquark 
limit by a = /?. 

A. Dipole Model 


We begin by examining a model dehned by the normalized momentum-space wave func¬ 
tion 


(j){abc) = 1 1 -h 

TT^ L J L 


-2 r 


-2 


(3,1) 


This form is motivated by its simplicity and by the fact that it guarantees that the result- 


Q. 


As we note in the following 


ing form factors fall off as 1/Q^, as predicted by QCD 
subsection, however, this requirement is not strictly required since we are working within a 
nonrelativistic framework that may not be applied reliably to large values of Q^. 

We have studied the integrals of Eq. (EH) numerically for arbitrary a/(3 and computed 
them exactly in the limit a//3—For particularly small diquarks, one hnds 

= (2 + COS20) • y , 


^Mp — 


5-1-3 cos 26 + sin 29 (3'^ 
2 -|- sin 29 2 

( 3 ^ 


^En = - COS 20 ■ y , 


^Mn — 


3-I- COS 20-I-sin 20 (3'^ 
' ' ~2 


(3,2) 


1 + sin 29 

In fact, one can shown that any ratio of the squared charge radii above is independent 
of the model wave function. Supposing that the 0 dependence of the charge radii varies 
weakly in a//?, and noting the basic pattern of the charge radii data, r% p ~ ~ r\p,^ 

and < 0, one may use Eqs. (El to predict disfavored regions of 0. In particular, 
r\n <0 only for | cos0| > 1/^/2':^ 0.7, while r'^^p/r'lpp& [0.75,1.0] for —0.4<cos0 < 1.0 and 
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FIG. 1: Plot of In(x^) for fixed mixing angle, assuming dipole wave functions. 


^ [0.93,1.33] for 0.4 < cos^ < 1.0. One therefore expects a favored region (noting 
that cos6' = ±l are the same solution) near | cos^l = 1. Indeed, this is borne out by our full 
numerical simulations. 

Numerical evaluation of the integrals in Eq. (EH) was performed using both Gaussian 
quadrature and adaptive Monte Carlo integration techniques. A fit was performed for 
the three ratios of the squared charge radii, rj^^p/r'^^p, ^E,n/^E,py '^M,n/'^E,pj while the 
overall scale was hxed afterwards using the experimental value of r‘^ p- (Nearly identical 
results were obtained when the overall scale was included in the function.) There are 
s.ve.al dete....aa«„as of fhe cha.go and .agnefic nadii in fhe Ufonatune H U £|. 
We have used the radii and uncertainties from the review of Hyde-Wright and de Jager 12j |. 
which are summarized in Tabled 

Values of a/P between 0.1 and 1.0 were sampled and the value of cos6' yielding the lowest 
was determined in each case. (For smaller a/f3, the asymptotic forms given in Eq. (j3.2jl 
yield a that is signihcantly larger than the results discussed below.) We hnd that the 
lowest x^ for C(/P G [0.1, 0.85] always corresponds to |cos6*| near one; this result is illustrated 
in both Fig. [Hand Table ITTl For any hxed choice of cos^ such that | cos6^1 > 0.85, one can 
hnd a value of a/(5 that provides a good ht to the data, i.e., the qualitative behavior of 
as a function of a/f3 is the same as the cos 6 = 1 curve in Fig. ^ Table iTTI illustrates that 
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TABLE I: Values for the nucleon charge and magnetic radii. 


Observable 

value ± uncertainty 

Reference 

rEp 

0.895 ± 0.018 fm 

M 

fMp 

0.855 ± 0.035 fm 

M 

' En 

-0.119 ±0.003 fm2 

[16] 

rMn 

0.87 ±0.01 fm 

m 


the lowest over a wide range of a//? (corresponding to some cos6' value near one) varies 
only weakly in a//3. The best fit point globally is at a//3 = 0.266 and cos^ = 0.917. In 
particular, for a = 0.244 fm and (3 = 0.917 fm, one hnds 

^Ep = 0.895 fm , VMp = 0.896 fm , 

VMn = 0.862 fm , r|„ = -0.119 W . (3.3) 


B. Gaussian Model 

The form of the momentum-space wave functions that we have just considered was mo¬ 
tivated by the asymptotic behavior of the form factors at large momentum transfer. How¬ 
ever, the charge radii are evaluated at zero momentum transfer and one might anticipate 
that other choices for the momentum-space wave functions (with differing high-Q^ behav¬ 
ior) should produce qualitatively similar results. We demonstrate this point by considering 
Gaussian wave functions. Although Gaussian wave functions lead naively to form factors 
that fall off too quickly at high this limit is also outside the range of validity of our 
nonrelativistic treatment; we do not let this point concern us further. We simply use the 
Gaussian case to illustrate the model-independence of our charge radii results. To proceed, 

a/(5 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 

cos6» 0.89 0.93 0.97 0.99 1.00 -1.00 -0.99 -0.88 
1.33 1.32 1.34 1.35 1.40 1.46 1.41 20.7 


TABLE II: Minimum of x^ in cos 9 for fixed a/f3, assuming dipole wave functions. 













FIG. 2: Plot of In(x^) for fixed mixing angle, assuming Gaussian wave functions. 


we take 



exp[-a^(pL“'“V] exp[-,3“(p!“'“>)"] . 


(3.4) 


Evaluation of the electric and magnetic charge radii proceeds exactly as described in Sec¬ 
tion m the only notable difference is that the necessary integration can be done more easily. 
In particular, the master integral of Eq. (EH) evaluates to 


J(7,5) = 


4a2/?2 


3/2 


X exp 


-1- cos^ 7 

[2(q;^-|- /3^) -1- (a^— /3^) (sin 27 sin 2 h— 2 cos ^7 cos 25)] 1 


(a^ -f 13‘^y — (a^ — 13‘^y cos^ 7 


J 


.(3.5) 


The qualitative features of the results are the same as in the previous case: Good hts are 
only obtained for cos 6 * near ±1. For example, for 0.2 < q;// 3 < 1, we find that the best fit 
values all occur for cos 6 > 0.85. This feature of our solutions is evident in Fig. |2l which 
shows In(x^) (using the data selection in Ref. |l^) as a function of a/f3, for three choices of 
cos 9. Only for cos 6 near one does the curve develop a deep local minimum that corresponds 
to a good fit to the data. For any cos^ > 0.85, there is an a/(3 that provides a low fit 
to the data. The set of these a/j3 values spans a wide range, as shown in Table HTTl below. 
As in our previous model, these best values remain relatively flat in a//? until one gets 
very close to a/(3 = 1 from below. Thus, our qualitative conclusion remains the same: the 
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a//3 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 

cos 9 0.85 0.90 0.95 0.98 1.00 1.00 -1.00 -0.99 
1.38 1.34 1.34 1.36 1.40 1.54 1.55 1.65 

TABLE III: Minimum of x^ in cos 6 for fixed a/f3, assuming Gaussian wave functions. 

current data does not strongly favor either the no-diquark or small-diquark solutions. The 
best fit point globally is in this case at a/f3 = 0.360 and cos6 = 0.927. For a = 0.311 fm 
and (3 = 0.865 fm, one finds 


^Ep = 0.895 fm , VMp = 0.897 fm , 

^Mn = 0.863 fm , = —0.119 fm^ , (3.6) 

where the overall scale of the parameters has been chosen so that reproduces the correct 
experimental central value. 

In both the dipole and Gaussian models we find that the most constraining feature of 
the experimental data is the smallness of the neutron electric charge radius relative to the 
other radii. The increase in x^ away from the best-fit solutions comes almost entirely from 
difficulty in accommodating r|;^. The value of this datum drives the fits to large values 
of cos 6, and toward particular values of a/13, where the model can best accommodate the 
two distinct length scales represented by the data. We hnd that our results do not change 
qualitatively if we change the experimental value of to another small number of either 
sign, or even set it identically to zero. 

It is also worthwhile to note that more distinct global minima can be obtained if one 
reduces the uncertainties on the charge radii. For example, if one arbitrarily shrinks the 
uncertainty of the ratio r\jp/r‘^p by a factor of two in our x^ function, one finds that a new 
global minimum appears near a//3 ^ 0.7, and that the depth of this minimum is enhanced 
in comparison to the results presented in Table IIIII (For example, the difference between 
the x^ at the minumum and the x^ at a//? ~ 0.2 is increased by an order of magnitude when 
the uncertainties are reduced in this way.) The significant spread in central values of the 
various existing estimates of the charge and magnetic radii ll^ llJ, ll^ suggests that an 
improvement in our input data is possible. In this case, a more definitive result on diquark 
correlations could plausibly be obtained from the charge radii alone, without consideration 
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of the higher momentum behavior of the form factors. 


IV. CONCLUSIONS 

The possibility of signihcant spatial correlations between quark pairs within three-quark 
baryons has been met with renewed interest in recent years. In this Letter, we have consid¬ 
ered whether precise measurements of the electric and magnetic charge radii are sufficient 
to diagnose whether any correlation identihable as a diquark is present within the nucleon. 
Since charge radii are evaluated at zero momentum transfer, we have considered nonrelativis- 
tic quark models that properly encode the total antisymmetry of the spin-flavor-momentum 
wave functions. Unlike some of the other diquark models in the literature, ours allows us 
to continuously interpolate between a small diquark (where such antisymmetrization is not 
necessary) and a mild spatial correlation that is not signihcantly smaller than the size of 
the nucleon. Working in momentum space, we have considered wave functions of dipole and 
Gaussian form. While the former has a more realistic high-Q^ fall-off, our results for the 
charge radii in each case agree qualitatively, suggesting that our conclusions are relatively 
model independent. We hnd a good £t to the data for a range of spatial correlations, pro¬ 
vided that the correlated quarks are primarily in a scalar diquark state, as we have defined 
it. Local minima of comparable are present for both small-diquark and no-diquark so¬ 
lutions, suggesting that a conservative assessment of the experimental uncertainty of the 
nucleon charge radii does not yet allow a definitive statement on the existence of a diquark 
in the nucleon. However, the spread in existing determinations of the charge radii suggests 
that there is room for a reduction in the uncertainty of our input parameters; in this case, 
it should be possible, in principle, to determine whether or not small diquark solutions are 
favored, without reference to the higher-Q^ behavior of the form factors. In the meantime, 
the overall results of our analysis suggest that a fully relativistic extension of our approach 
{i.e., evaluation of the nucleon form factors assuming properly symmetrized wave functions 
that can smoothly interpolate between spatially correlated and uncorrelated quark pairs) is 
well motivated and worthy of further investigation. 
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